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Analog PLL's are being replaced by a pure digital means of generating any desired frequency. 
Computer circuits can generate millions of samples of sine waves which are converted to an analog 
wave form using a digital to analog converter. The process is almost the reverse of how a CD player 
works. The number of samples per second required is determined by the Nyquist theorem. Nyquest 
states that to re-construct a given wave form at least two samples for each 360 degrees of the wave form 
are required. Simply stated, the sample clock must be twice the highest desired frequency of the 
desired wave form. The device that performs this magic is called a Direct Digital Synthesizer or DDS. 
The DDS consists of a digital accumulator, a wave form fiinction generator, and a digital to analog 
converter. The digital accumulator is an adder and latch circuit. On each sample clock the accumulator 
adds the previous sample to a constant. The value of the constant determines the frequency of the wave 
form generated. For example, if we used a 16 bit accumulator, and the constant was 256, it would take 
256 clock samples to complete a full cycle. I.E.: 65536 / 256 = 256. The output frequency would 
therefore be sample clock freq / 256. The wave form function generator takes the output value of the 
digital accumulator and transforms it into a value of the wave form at a given point in time. In most 
cases the wave form function generator is a read only memory containing a sine function look up table. 
The sine value of the digital accumulator is then presented to the digital to analog converter. It is now 
only necessary to remove alias products from the wave form by passing it though a low pass filter. 
Several corripanies now make DDS components on a single chip. Analog Devices has several DDS 
devices complete with computer interfaces, clock generator, and digital to analog fiinctions on a single 
chip that will function to as high as 300mhz. 

Special microprocessors (digital signal processors or DSP) are now fast enough to analyze and modify 
wave forms in real time. If an analog signal is first passed though a low pass filter to bandwidth limit 
it, and then to an analog to digital converter such a microprocessor may operate on the resulting data 
and then output the modified data to a digital to analog converter. After again passing though a low 
pass filter the resulting signal is re-created, but in modified form. By the use of proper algorithms, it is \ 
possible to filter out QRM and QRN (interference and noise) from the signal. Several filters may be 
applied at once to the signal. If the DSP and analog converters are fast enough, and a low frequency IF 
is used, the IF signal may be processed and the signal can be demodulated in the digital domain. 

In the early days of amateur single side band two methods of signal generation and demodulation were 
considered. The filter method used a high frequency filter (crystal or mechanical) to remove one of the 
unwanted side bands and to help suppress the carrier. The phasing method used a special mixer circuit 
(balanced modulator or image reject mixer) to remove the unwanted signal components. The filter 
method won out and is used today. However it is now possible to make use of the phasing method with 
digital means. The phasing circuit required that the audio signal and rf drive into the mixer (or doubly 
balanced modulator) be split into two 90 degrees apart. A DDS circuit with both sine and cosine 
function generators will generate precisely split quadature signals. The audio signal may be precisely 
split by use of a DSP. The doubly balanced modulator will work in reverse to demodulate a SSB signal 
and to select the desired side band. In this case the circuit generates a quadature audio signal pair, and 
the DSP must combine them into a single signal for only the desired side band. It is now possible to 
build a direct conversion SSB set which will generate and receive signals directly at the air frequency. 
Such a set has been described in QST, but it used analog components to generate and combine the 
quadature signals. Who will be the first to build such a digital radio? 

As DSP circuits improve, and get faster a completely software radio is possible. Such a rig would have 
the antenna connected (via a tuned front end and perhaps some pre amplification) to an analog to digital 
converter. Tuning and demodulation would be done completely in software, and the output of the 
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digital to analog converter presented (after a low pass filter) to the speaker. While such a radio is still 
in the future, I suspect it will find it's place in the ham shack of many of today's amateurs. 
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1 Introduction 

This project will explore issues associated with implementing a OMSK (Gaussian Minimum 
Shift Keying) transmitter that is intended for use in GSM cell phones. This communication 
standard encodes information using phase modulation of a constant envelope sine wave 
signal. By constraining the transmitter output to have constant envelope, a nonlinear power 
amplifier can be used that broadcasts the signal with high efficiency (usually around 50% 
efficiency, as compared to about 10% efficiency encountered with linear power amps). 

A proposed transmitter for this task is shown in Figure 1, which employs a PLL that 
contains an I/Q modulator in its feedback loop. As you proceed through the project, you 
will (hopefully) come to understand the value of placing the I/Q modulator in the PLL 
feedback loop rather than directly using its output as the transmitter output, and develop 
a good understanding of RF transceivers in general. 

Figure 2 displays the proposed receiver used to select and demodulate the desired GMSK 
RF signal, which is implemented using a direct conversion architecture. For simphcity, we 
will ignore the many difficulties encountered with this structure, such as high sensitivity 
to DC offsets, antennae impedance variations, and local oscillator feedthrough. The noise 
source in the figure represents the overall input-referred noise of the receiver (it is not, of 
course, purposefully added!). 

2 Background 

Two common phase/frequency modulation schemes are in use today — GMSK and GFSK 
modulation. GMSK modulation focuses on phase as the modulation variable, while GFSK 
(Gaussian Frequency Shift Keying) focuses on frequency. It turns out that GMSK is superior 
in terms of its spectral efficiency, but GFSK is easier to implement. As a result, GMSK is 
typically used for cell phones (using the GSM standard), while GFSK is often employed in 
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Figure 1: A GMSK modulator implemented as an offset PLL. 
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Figure 2: A direct conversion receiver for GMSK detection. 

cordless phones (using the DECT standard). Note that cordless phones can get by with much 
lower performance (and lower cost) than their cell phone counterparts for reasons discussed 
in class. 

Both modulation methods use a Gaussian transmit filter to smooth the data signal, as 
illustrated in Figure 1, in order to reduce its spectral content. This filter is implemented 
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in discrete-time with a sample period of T, and will be denoted as P[k]. If the input data 
pattern consists of modulated, unit area impulses with symbol period (T^) spacing between 
them, then the transmit filter is formulated as P[k] = P{Tk) such that 



P{t) = h 



T 1 
2 y/2^a 




*Tect{Td,t) 



(1) 



where denotes convolution, and 



a = 



.833Trf 



rec 



±{Td,t) = 



0, elsewhere. 



In the case where the data is "square- wave" in nature, as shown in Figure 1, then one 
simply removes the rect function (and the associated convolution operation) in the above 
formulation (note that you will need to appropriately scale the response, as discussed later). 
In either case, P[k] is parameterized by two characteristics: 

• BTa: ratio of bandwidth of P{e^'^''^^) to data rate 

• h\ modulation index, ^defined as: 



For the popular standards of GSM and DECT, we are required to have: 

• GSM: BTd - 0.3, h = 0.5 

• DECT: BTd - 0.5, /i = 0.5 ± .05 

Therefore, both specify the modulation index to be h=0.5 — for GSM, this value must be 
exactly maintained, while for DECT it can vary ±10% about that nominal value. However, 
GSM requires the BT^ product to be much lower, which implies a much better spectral 
efficiency (i.e., required bandwidth for a given data rate). 

The impact of the diflFerent BTd products for GSM and DECT can be viewed in the 
time or frequency domains. Here we focus on the time domain, and use eye diagrams to 
convey the key issues. Figure 3 illustrates ideal eye diagrams for the instantaneous frequency 
and phase signals corresponding to a GSM transmitter with BTd — 0.3. In this case, we 
are interested only in phase (since GMSK is used), and therefore the I/Q diagrams are of 
primary interest. Note that the eye for instantaneous frequency is relatively poor compared 
to its I/Q counterpart — extraction of the modulation information based on frequency would 
therefore yield rather poor performance (which is further exacerbated by issues related to 
synchronous vs asynchronous demodulation, a topic we will not touch on here). 

Figure 3 illustrates ideal eye diagrams corresponding to a DECT transmitter implement- 
ing GFSK with BTd = 0.5. For GFSK, we are concerned with the instantaneous frequency 
rather than the I/Q components. Note that the eyes for instantaneous frequency (as well as 
its I/Q counterpart) are quite open compared to the GSM modulator. We therefore see that 
high vales of BTd lower the level of intersymbol interference at the cost of poorer spectral 
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(a) (b) 

Figure 3: Eye diagrams correponding to GMSK modulation with BTd = 0,3: (a) eye diagram 
of instantaneous frequency, (b) eye diagrams of I and Q components 

efficiency. A final issue is that one would observe the I/Q eye diagrams undergo catastrophic 
degradation should the value of h vary even slightly from 0.5, while the instantaneous fre- 
quency eye diagram would simply change in amplitude while still remaining wide open. Since 
frequency is the signal of interest in this case, GFSK systems allow such sloppiness in the 
control of h, whereas GMSK system are intolerant of it. 




(aj (b) 

Figure 4: Eye diagrams correponding to GFSK modulation with BTd ~ 0.5: (a) eye diagram 
of instantaneous frequency, (b) eye diagrams of I and Q components 

The transmitter structure shown in Figure 1 achieves accurate phase modulation, and 
can be used for either GMSK or GFSK systems. However, the abihty of GFSK systems to 
endure sloppiness in their modulation index, /i, allows for simpler transmitter structures. In 
fact, it is not uncommon for GFSK modulators to be implemented by directly modulating 
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aVCO. 

While the GMSK receiver shown in Figure 2 appears to be quite simple, it is quite chal- 
lenging to implement in practice. In comparison, Figure 5 illustrates a GFSK receiver based 
on a phase-locked loop architecture. Here the receive VCO tracks the phase of the input 
signal after it has been filtered and passed into a limit amplifier, and the resulting instanta- 
neous frequency signal is obtained by observing the input voltage to the VCO. Though this 
design may appear to be more complex than its GMSK counterpart, it is easier to implement 
in practice. 
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Figure 5: A GFSK receiver based on a PLL frequency descriminator. 



3 Tasks 

We now turn our attention entirely to GMSK, and focus first on the design of the transmitter 
structure shown in Figure 1. We will do so in steps, examining first the I/Q generation, then 
the PLL design, and then simulation of the overall system. We will then turn our attention 
on the receiver, particularly with respect to its simulation. In all cases, you should assume 
the following: 

• Modulation data rate (l/T^) = 1 Mbit/s 

• Reference clock rate (l/T) - 100 MHz 

• Center frequency of RF output spectrum fjip — 900 MHz 

• Tristate PFD is used for the PLL, 

• VCO has no parastic poles, so that its small signal voltage-to-phase model is simply 

• VCO noise (output referred) has a slope of -20 dB/dec over the frequency range of 
interest, with a value of -140 dBc/Hz at 10 MHz offset 

• Simulation sample rate — 9 GHz. 

Note that, in practice, the GSM standard actually has a data rate of roughly 270 kbit/s — 
we are striving to increase it by a factor of 4. 
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1. In this task we will implement the I/Q generator block shown in Figure 6 using the 
CppSim simulator. The following subsections will guide you through this process. 
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Figure 6; digital I/Q generator. 



(a) Implement the data generator and Gaussian transmit filter in CppSim, and plot 
the resulting eye diagrams produced by a data sequence consisting of 30 symbols. 

i. Use Matlab to compute the filter tap coefficients of the Gaussian transmit 
filter over 5 symbol periods. Assume the sampling period of the tap coeffi- 
cients corresponds to that of the reference frequency, and that the input data 
sequence wifi be "square" wave in nature (i.e., don't convolve with the rect- 
angle function, but you do need to figure out the appropriate scaling factor 
based on the reference frequency and symbol rate value. To do so, recall that 
the peak-to-peak frequency deviation that occurs after passing the prbs data 
sequence into the Gaussian filter must match that specified by h and T^. As- 
sume that the prbs data sequence alternates between 1 and -1 and is sampled 
according to the reference frequency). Plot the filter impulse response over 5 
symbol periods using the stem command in Matlab. 

ii. Create the prbs data source module in CppSim using the Rand class ("bernouUi" 
behavior). To do so, you will need to clock the Rand object according to an 
input clock (which will be fed by the reference frequency) such that it updates 
its output value at the correct data rate. Note: see the next part to see how 
to "clock" objects. Plot the output of the module and its eye diagram over 30 
symbols (you can use the vco module to create the 100 MHz reference clock). 
Note that the output should alternate between 1 and -1. 

iii. Create the Gaussian transmit filter module in CppSim. To do so, save the 
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filter coefficients computed in Matlab to a file (make sure that you save in 
ASCII format, and that the coefficients are in "vertical" format in the file 
(i.e., one coeflficient per line)). To give you an idea of how to realize the 
module code in modules. par, consider the following example: 



module: sqrt_raised_cos-filter 
parameters: 

inputs: double in, double elk 
outputs: double out 

classes: Filter mt("l"."l"). List listl(), EdgeDetect clkedge() 

static, variables: 

init: 

list 1 . load ( " r CO s f 1 1 . dat " ) ; 
filt.set(listl,"l"); 

code: 10 
if (clkedge.inp(clk)) 



Plot the output of the cascaded data source and Guassian filter, and its eye 
diagram, over 30 symbols. Label the peak-to-peak frequency deviation and 
symbol period in the eye diagram. 

(b) Implement the remaining I/Q generator blocks, simulate its output over 30 sym- 
bols, and then create an overall symbol for it. 

i. Compute the value of Kp^ in the frequency-to-phase converter block. To do 
so, assume that the sample period corresponds to the reference frequency, 
1/T, and recall that z corresponds to e^^'^ and phase is related to frequency 
by the equation 



ii. Implement the frequency-to-phase converter block and the cosine and sine 
converter functions in CppSim (combine the cosine and sine converters in one 
module). Again assume that all blocks are clocked by the reference frequency, 
and be sure to use the Filter class to implement the accumulator function for 
the frequency-to-phase converter block. Plot the I and Q signals, along with 
their eye diagrams, over 30 symbols. 

iii. You have now completed the digital portion of the I/Q generator, so place 
each of your blocks (though not the reference frequency source) into a schematic 
called dig_iq-generator that has a clock input and I/Q outputs. For the ex- 
ercises that follow, note that the D/A converters do not need to be added in 
CppSim since our digital blocks are implemented in analog fashion, and the 
zero-order hold function is implemented intrinsicly. 

(c) Now put RC lowpass filters at the outputs of the I and Q outputs of the I/Q 
Generator (note that these will *not* be clocked by the reference). Set the cutoff 



{ 

filt.inp(in); 
out = filt.out; 



} 
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frequency of the filters first to 500 kHz and then to 10 MHz, and plot the resulting 
I/Q diagrams at the filter outputs for each case. Why does the 500 kHz version 
look distorted? Assume that the cutoff frequency will be 10 MHz for the following 
task sections. 

2. We will now focus on modeling and design of the PLL portion of the transmitter shown 
in Figure 7. You will need to use the PLL Design Assistant program, as used in HW6, 
to complete the design exercises. 




Figure 7: PLL portion of the GMSK transmitter. 



(a) Model the relationship between the phase signal, the phase of the VCO 

output, ^out{t)i and the output phase of the I/Q modulator block, ^out2{t), as 
labeled in Figure 8. Do this in two steps. First, express the VCO output as a sine 
wave, and compute the resulting I/Q modulator output in terms of a sine wave. 
Based on this result, draw a block diagram of the relationship between the phase 
signals. 




Figure 8: Analog I/Q modulator block. 
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(b) Based on your block diagram in the previous part, draw the Unearized, frequency- 
domain block diagram model of the PLL system shown in Figure 7. Provide an 
expression for the open loop gain function, A{s). For parts to follow, we will 
define the closed loop function G{f) to be 

G(/)= 



(c) Compute the transfer function from the phase signal, to the VCO output 
phase, ^out{t) in terms of G{f). We will refer to this function as B{s) in the parts 
to follow. 

i. Is B{s) lowpass, bandpass, or highpass? 

ii. If we instead look at the sine wave signals associated with out{t), etc., rather 
than their respective phases, is the filtering action performed by the PLL 
lowpass, bandpass, or highpass? What are the advantages of using the PLL 
to implement the filter over simply using a passive filter? 

(d) Given your answer in the previous part, do you see any advantage of incorporating 
the I/Q modulator in the feedback loop of the PLL rather than simply using the 
I/Q modulator as the transmitter (i.e., with a PLL only being used to provide a 
static local oscillator signal)? Consider noise, as well as mixer nonidealities such 
as LO feedthrough, in formulating your answer. 

(e) Design the PLL such that the B{s) transfer function corresponds to a second 
order Butterworth filter response with a bandwidth of 2 MHz. Assume a type I 
implementation for the PLL, a charge pump current of 100 /xA, and a Ky value 
of 30 MHz/V. In particular, compute the divide value, N, the loop filter transfer 
function, H{s), and the PFD gain factor, a. 

(f) Repeat the previous design exercise for the case where B{s) corresponds to a 
second order Butterworth filter response with a bandwidth of 4 MHz that is 
cascaded with two real-valued parasitic poles placed at 12 and 17 MHz. 

(g) Calculate the input-referred voltage noise to the VCO that yields to its phase 
noise conforming to the specifications given at the beginning of the Tasks section 
(i.e., -140 dBc/Hz at 10 MHz offset). - 

3. We will now implement the PLL circuit in CppSim and combine it with the I/Q 
generator module. This task should be fairly straightforward given your experience in 
Homework 6, but there are a few hurdles to overcome that we will step you through. 

(a) Implement the I/Q modulator shown in Figure 8 in CppSim by creating ideal 
mixer blocks and the hmit amp module, and then combining those modules with 
the adder module in a schematic/symbol called phase_shifter. Assume that the 
I/Q modulator has two inputs corresonding to the VCO output and its quadrature 
component (i.e., don't implement the 90 degree phase shifter). 
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i. When creating the hmit amp module, you must understand that the output of 
the I/Q modulator is essentially a sine wave, but that the input to the divider 
assumes a square-wave with interpolated edges as discussed in the paper 
located at /mit/6.976/CppSim/Doc/paper.pdf. The following is suggested 
code to implement the hmit amp — explain what the code is doing. 

module: phase_shifter_ limit _amp 

parameters: 

inputs: double in 

outputs: double out 

classes: 

static- variables: double prev.in, int state 
init: 

state = -2; 
prev.in = 0.0; 

out = -1.0; 10 
code: 

if (state == -2 in >= 0.0) /* positive transition */ 

{ 

out = (in+prev_in)/(in— prev-in); 
state = 1; 

} 

else if (state == 1) 
{ 

out = 1.0; 

state = 2; 20 

} 

else if (state == 2 iiii in < 0.0) /* negative transition */ 

{ 

out = — (in+prev_in)/(in— prev-in); 
state = — 1; 

} 

else if (state == -1) 
{ 

out = -1.0; 

state = -2; 30 

} 

prev-in = in; 



(b) Implement the VCO by modifying the vco module already included in the Cpp- 
Modules library such that it now has quadrature sine wave outputs (i.e., cosine 
and sine functions of the VCO phase). These outputs will connect up to the 
inputs of the I/Q modulator just implemented. 

(c) Now implement the overall PLL in CppSim, and choose the various parameters so 
that the dynamics conform to those calculated in Task 2, part (e) (do not include 
noise sources). Include the I/Q generator from Task 1, and name the overall 
schematic plLtransmitter. Simulate the overall system over 30 symbols and plot 
the voltage at the input to the VCO and its eye diagram (you may want to cut 
out the initial transient when plotting the eye diagram). Relate this signal to the 
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instantaneous frequency waveform that you should expect to see. Troubleshoot 
your simulation if it does not agree with what you expect to see. 

(d) In Cadence, copy plLtransmitter to a new schematic called plLtransmitter2. Ad- 
just the parameters in plLtransmitter2 so that its dynamics conform to those 
calculated in Task 2, part (f). As in the previous part, simulate the overall sys- 
tem over 30 symbols and plot the voltage at the input to the VCO and its eye 
diagram. Do you see any change in performance in the eye diagram (i.e., is it 
better or worse?)? Why did it change? 

(e) Plot the output spectrum of the plLtransmitter module using the comp.psd Mat- 
lab script as done in Homework 6 (note that you may need to adjust some things 
in the script to match the given system). You will want to dramatically increase 
the number of symbols to achieve a nice plot. Also, you will need to take advan- 
tage of CppSim's abihty to decimate signals before they are saved in the output 
file in order to simulate over a long enough time span — see the description of 
the "output:" command on pages 32-33 of the CppSim manual. 

(f) Add VCO noise according to your calculations in Task 2, part (g), and then 
plot the resulting output spectrum as done in the previous paxt. Compare the 
spectrum before and after adding VCO noise — do the results make sense? Note 
that, in practice, there would also be detector noise in the PLL and quantization 
noise from the D/A. We have ignored such sources in the current assignment for 
the sake of simplicity, but they would be fairly straightforward to add to the 
simulation. 

4. We now examine the receiver, and will implement a simplified version of it shown in 
. Figure 9 in the CppSim framework. 
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Figure 9: Simplified GMSK receiver to be implemented in CppSim. 



(a) Implement the receiver shown in Figure 9. Assume that the noise source has zero 
variance (for now), and that each channel select filter is simply a pair of cascaded 
first order filters with a cutoff frequency of 5 MHz (DC gain of 1). Create an 
overall Cadence symbol for the receiver called simple Jq_receiver that has 3 inputs 
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(to receive transmitter output and the quadrature sine waves at frequency fnf) 
and two outputs (the baseband I/Q signals). 

(b) Now create a schematic that combines the plLtransmitter and simple_iq_receiver 
blocks and also supplies the required LO frequency to the receiver. We will assume 
an ideal channel, so simply connect the transmitter output to the receiver input. 
To create the quadrature LO signal, use the VCO block you created in Task 3 that 
produces quadrature sine waves. Simulate the system over 30 symbol periods and 
plot the resulting I/Q signals at the receiver output along with their eye diagrams. 
Note that you may need to adjust the phase of the LO signal to the receiver in 
order to get your eye diagrams to look right — explain why that is the case (hint: 
think about the I/Q signals in terms of constellation points). 

(c) Compute the value of the receiver noise variance that will yield an SNR of 10 dB 
at the I/Q outputs of the receiver (assume that the noise is" added as shown in 
Figure 9). Here we define SNR as the ratio of the spectral density of the received 
signal to the spectral density of the noise at low frequencies, as shown in the 
figure. Simulate the overall system with the included noise and plot the resulting 
I/Q eye diagrams. 

(d) Given the added noise of the receiver, what are the tradeoffs between having a 
large versus small bandwidth for the receiver channel select filter? 

5. We have so far assumed that we only want to modulate in a channel that is centered at 
900 MHz. Suppose that you wanted to use the same RF transmit architecture (with a 
constant reference freqeuncy of 100 MHz), but wanted to be able to change the center 
frequency in 1 MHz steps. Explain how you would modify the transmitter architecture 
to accomplish this task, and prove it by simulating your modified structure in CppSim 
with fRF = 901 MHz. 
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